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Abstract. Every isotropic probability measure on the space R°° of real sequences 
x=(xi,x 2i ...) is a convex combination of the measure concentrated at 0 and a member of / 0 (/?“), 
the set of all isotropic probability measures p„ on /?°° with p„({0})=0. Each p„e / 0 (/?°°) is 
completely determined by any one of its finite-dimensional marginal distributions p n . Each p n 
has a density function /„ with dp n (x x x n )=dx x • • • dx n f n (x? + ■■• +x„ 2 ). Each f n is com- 

pletely monotone in 0<£ <«» (hence analytic in the right complex £ half-plane), and 

oo 

JC nn r(n/2y l J d£Z n,2 - l f n (£) = 1 . 

o 

Every / which satisfies these two conditions is /„ for a unique p„e /(//?“)• Hence the equation 

oo oo 

n jdCf 2(D = J«)e-' ? 

defines a bijection between 1 0 (R°°) and the set of all probability measures p on 0 < r < °°. If 

oo 

p„e / 0 (R “) then pJi{x : X x,- 2 <«>})=0, so p x is not a "softened” or "fuzzy" version of the ine- 

1=1 

oo 

quality X x,- 2 < 1. If the prior infomiation in a linear inverse problem consists of this inequality 
1=1 

and nothing else, stochastic inversion and Bayesian inference are both unsuitable inversion tech- 
niques. 
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Introduction. Let R be the real numbers, R n the linear space of all real n -tuples, and R°° 
the linear space of all infinite real sequences x = (x u x 2 ,.. )- Let P n :/?*"->/?" be the projection 
operator with P H (x) = (x x ,...,x n ). Let p„ be a probability measure on the smallest <j-ring of sub- 
sets of R” which includes all of the cylinder sets P„ l (.B n ), where B n is an aibitrary Borel subset 
of R n . Let p n be the marginal distribution of p„ on R H , so p n (B n )= p„(P~ l (B „ )) for each B n . 
A measure on R n is "isotropic" if it is invariant under all orthogonal transformations of R * . The 
measure p„ will be called isotropic if all its marginal distributions p n are isotropic. The set of all 
isotropic probability distributions on R” will be written /(/?“). The present note describes all 
members of /(/?“)• The result calls into question both stochastic inversion and Bayesian infer- 
ence, as currently used in many geophysical inverse problems. 

Necessary Conditions for Isotropy. Let 0=(0,0, ••• ) and let p2, be the member of 
/(F~) such that ({0})= 1. If p„e anda,/?£0 anda+/J = l, then ap..+/Jp ° e I(R°°). 
Conversely, if and /><«,( {0})=/J, then p, 0 =(\-fi)p 0C +Pp® where p„s I (R°°) and 

/U{0})=0. Therefore it is necessary to study only those p^e I(R°°) for which p oo ({0})=0. 
They constitute the subset / 0 (F °°) of I (R °°). 

If p„ e I 0 (R “), for every £ in 0 <£ < <» define 


F n (£)=pM*'-Xi + ••• +x?>£}). [1] 

Then F n is right semi-continuous, and 

^„(0)= 1 [2a] 

F„(oo)=limF„(g) = 0. [2b] 

Also, if n <N and a<A , then 

O^F^A^F,,^)^^!. [2c] 


Properties sufficient to characterize the members of 7 0 (/?°°) are given in 
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Theorem 1: Suppose p«, e / 0 (^°°) and F n given by [1]. Then for each integer n £ 1, F n (£) is ana- 
lytic in the open right half plane of complex There is a function /„(<!;), also analytic there. 


such that for every Borel subset B n of R " 

Pn(B„)= J dx x • dx n /„(x 2 + •■■ +x n 2 ). [3a] 

B. 

In particular, if 0 <a < <*> then 

oo 

F n (a) =n na T(nl2T x J d^ n!2 ~ l f n (§) . [3b] 

a 

The /„ arc related by 

oo 

/„(0 = |rfr7(n-^rV» + i(n) [3c] 

oo 

| dT)(n-$r' A f H (77) [3d] 

oo 

/„(£)=* Jd 77 / B+ 2 ( 77 ) [3e] 

/» +2 <5)=-* _ V»<£> t 3f i 

For every /? in 0 </? < » 

lim (9) = 1 . [3g] 

n — >«» 


PROOF: Let S(n-l) denote the unit sphere in R n , and let IS(n-l)l be its (n-l)-dimensional 
Euclidean content, 7:t nl2 T{n /2) -1 . Let IS(n-l)l 0 „(w) be the content of the part of S (n-\) 
where x„ 2 < 1 -w. Then 

w 

^ + l(w) = l- I5(n-l)l ISCn )!- 1 JrfCC" / 2 (K)' W - 

0 

Since p„ is the marginal distribution on R n of p „ +1 on R n+1 , 

oo 

F b (O = -|^„ + i(77)0„ + i(^). [4a] 

the right side being a Stieltjes integral. For any/? and £ satisfying £</?<£, d n <p n+ \(^h]) is con- 
tinuous in/? <77 <B, so integration by parts (1) permits the conclusion 
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j dF n+x (j))<l> n + x (gh]) + j diiF n+l (T 1 )d r fp n+x (l;fn) 

= Fn+i(B)<t> n+x (g/B)-F H+x (p)<l> n+ i<g/p). 

Here let/3 -»£+ and B The integrated parts tend to zero, so the Lebesque bounded conver- 
gence theorem permits [4a] to be rewritten 


|-” /2 f n (0= is(n-i)i \s(n)\~' ^dnn^+WF^Tixn^r*. 
Iterating this formula once, reversing orders of integration, and invoking the identity 


jd/KC-r?) 


leads to 


S-* a F n <§) = (.mi) | dCC <n+m F n+2 (D • 

By induction on n , it follows that F n (£) is infinitely differentiable in 0 <£ < °°. If we define 


[4b] 


f n (0=-*~ n,2 r(n/2)Z 1 - n %F n (g), [5a] 

then /„ is also infinitely differentiable in 0<<f; <<*> and [2b] yields [3b]. Then [3a] follows by 

straightforward integration theory. Then the definition of marginal distributions implies 


fn(X\+ • • • +X„) = J dx n+x f n+x (xf+ - +X„ 2 +1 ) , 

-oo 

which is [3c] with£ =x x + ■ ■ ■ +x n 2 ,T]=x x + • • • +x? +x . Also, 

oo oo 

fn(*\ + ■ +X„) = J dx„+ j | dx n+2 f n+2 (x l + " +^+2>. 


[5b] 


[5c] 


which is [3e]. Then [3f] follows from [3e], and [3d] follows from [3f] and [3c] with n replaced 
by n— 1. To prove analyticity, note that if q is an integer >0 and if 0<« </3, then by Taylor’s 
theorem with remainder 

F 2 «X) - F 2 <P) = £ (~d^F 2 (p) + -jy J d^-a)H-di) q+l F 2 (0 . [6a] 

i=l 1 • Q • a 
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But (-d^Y F 2 -n‘ f 2 i . so by [3b] 
t P 

~ J dSP(-%V* lF &) = F 2q+z( a ) - ^2, +20) ■ [6b] 

Hence, the Lebesque bounded convergence theorem implies that as a — >0 the integral in [6a] 
converges to 1 -F 2^+20)- Therefore 

(-^)‘>2<P) • [6c] 

i - 1 1 * 

All terms in the sum [6c] are nonnegative, and + 2 (jS) < 1 , so the series 

5 F ^' ><iS) t6d] 

converges absolutely (here F ^ =d£F 2 ). Therefore, the power series for F 2 (£) at| =/3 converges 
absolutely for all complex £ in the closed disk l£-/3 1 </J. Since p is arbitrary, F 2 (i;) is analytic 
for all complex with positive real part. By [5a], so is/ 2 (£) and then by [3c,d] so is /„(<!;) for 
every n >1. Hence so is F n (g) for every «>1. Furthermore, since [6d] converges, Abel’s 
theorem (2) implies that 

F 2 (0) - F 2 (fi) (^)'>2<J3) • [6e] 

«=i i! 

Together, [6e], [6c] and [2a] imply [3g]. 

COROLLARY 1: If one of the marginal distributions p n is known, is completely determined. 

00 

COROLLARY 2: Let H(a) be the set of x in R°° with X jc, 2 <a. Then p„(H ( <~)) = 0. This fol- 

1*1 

lows immediately from [3g] and the fact that H (°°) is the monotone limit of the sets H (a) (3). 

Sufficient Conditions for Isotropy. Let M(n) be the set of infinitely differentiable real- 
valued functions/ on the open half-line 0<^ <<» such that 

00 

K na x \n a r 1 J d^ nl2 ~'f <£) = 1 [7a] 

0 
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and also for every integer < 7^0 and every £ in 0 <£ < «> 

0 . [ 7 b] 

Note that if p„e I oC/?°°) and f n comes from p„ via [3a] then/„ e M(n). The converse is also 

true, and to prove it we need 

LEMMA 1: Suppose n £ 1 and/ e M(n). Then 


lim £ b/2 /<£) = 0 

[ 8 a] 

lim S" 72 / <£) = <) 
f-*o 

[ 8 b] 

00 

f(g) = ^dTl[-d t ,f(ll)] 

[ 8 c] 

OO OO 

(nil) j d^ n,2 ~ i f (§) = j d^ nl2 [-d^f (0) 
0 0 

[ 8 d] 

e M(n+ 2) . 

[ 8 e] 


PROOF: Let m=nl 2-1 and let 0 <a < A < ~. Integration by parts gives 

A A 

{m+\)\d$$ m f(£) = A m+1 / (A ) -a m+} f (a) + jd£ € m+1 [^/ <£)] . [9a] 

a _ a 

Fix a. The integral on the right in [9a] increases asA—>°° and yet is bounded, so it has a limit. 
Therefore lim A m+l f(A ) exists. By [7a] it cannot be positive, so we have [ 8 a], and hence [ 8 c], 

A — **> 

and also 

00 00 

(m+ 1 ) J dZS m f®=-a m+i f(cc) + Jd« m+ 1 [-V€)]. [9b] 

a a 

As a decreases to 0, the integral on the right in [9b] increases, and that on the left has a finite 
limit, so a m+l f (pc) approaches either +«> or a nonnegative limit. Then [7a] requires [ 8 b], and 
[9b] converges to [ 8 d]. Then [ 8 e] follows from [ 8 d] and [7b]. 

Now we can prove 

THEOREM 2: Suppose n is a nonnegative integer and / e M(n). Then there is a p„e 1 0 (R°°) 
whose marginal distribution p n on R n is given by [3a] with f„ =f . 
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PROOF: For every integer q £ 0, define /„ + 2 q (£ ) =7T _<? (-d{ ) q f(£)- liN-n is a nonnegative even 
integer, induction on [8c] implies 

oo oo 

/n(*i 2 + • • • +*w) = J <&/v+i J dxN+if n+2( x i + - • * +*n+ 2) • [10a] 

If N-n is a nonnegative odd integer, define f N from f N+l via [3c]. Then 

oo 

/n( x I + ' ■ ' +X N ) = J dXN+lf N+l( x l + ' ■ ' +*N+1 )• [10b] 

That [10b] also holds when N-n is nonnegative and even follows from [10a]. Therefore [10b] 
holds for all N ^.n. Use it inductively to define f N for 1 £ N < n . For N =n, [7a] implies 

j dx { ■ dx N f N (xf + ■ ■ ■ +X[})=1, [10c] 

and then [10b] implies [10c] for all N £ 1. Thus the probability distributions p N on R N given by 
f N via [3a] satisfy the Kolmogorov consistency condition. Then the existence of p„ follows 
from Kolmogorov’s Fundamental Theorem (4). 

COROLLARY l:\if e M(n ),/(£) is analytic in the open right half-plane of complex £ . 

oo 

COROLLARY 2: The equation F 2 (^) = ^dp(t)e~^‘ furnishes a bijection between the members of 
I o (R °°) and the probability measures p on 0 < t < «>. 

PROOF: Demanding that / 2 e M (2) is equivalent to demanding thatF 2 (^) be completely mono- 
tonic on 0<Z <°° (5). 

Examples and Applications. Setting f $>)=n~ x e~^ gives /„ (£ )=K~ na e^. This p M is the 
gaussian with independent x it x 2t ~, each having mean 0 and variance 1. Setting 
y 2 (^)=^ _1 v[^ v_1 -(l-^) v_1 ] with 0<v < 1 gives a p„ for which lim /„(£)=<» if n <2 and also if 

n = 1 and Vi<v < 1. Thus the densities /„(£) need not remain finite as£ — >0. 

The geophysical application is to inverse theory. An infinite dimensional linear space X of 
earth models x is given, along with a finite number of linear functionals, gj.X—>R, 


September 4, 1987 



George Backus 


9 


Proc. Natl. Acad. Sci. USA (1987) 


j=l,...,D+l. An observer measures D datay,- =g i (x E )+e i for i=l,...,D. Here x E is the correct 
earth model and e,- is the error in observing y,-. The observer wants to predict the value of 
z =g D+ i(x £ ). Since dimX =«, the problem is hopeless unless g D+l is a linear combination of 
SU‘—8d • or unless the observer has some prior information about x E not included among the 
data (6,7). One common sort of prior information is a quadratic bound on x E , a quadratic form Q 
on X such that x E is known to satisfy 

Q(x e ,x e )< 1. [11] 

Often [1 1] is a bound on energy content or dissipation rate (8). In stochastic inversion and Baye- 
sian inference, such a bound is often "softened" to a prior personal probability distribution p„ on 
X (8-10). In practice, X is truncated to an R n , and p n is used in the inversion. 

To see why this process is questionable, complete X to a Hilbert space with the inner pro- 

oo 

duct x x'=Q(x,x'). Let £i,£ 2 > ' be an orthonormal basis for X, and write x= £ x t £,• . 

i =1 

Then X becomes the subset H (°°) of R°° defined in corollary 2 to theorem 1. The prior informa- 
tion [11] can now be written 

ixfzi. [ 12 ] 

< =1 

If the observer wants to soften [12] to a probability distribution p w without introducing new 
information not implied by [12], then clearly he should take p„e I (R°°). He is unlikely to assign 
nonzero probability to 0, sop 1 But then pJX )— 0 by corollary 2 to theorem 1. Any 
prior personal probability distribution obtained by softening [12] without adding new information 
must deny [12] with probability 1. 
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